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ABSTRACT 

With the rapid progress in metallicity gradient studies at high-redshift, it is imperative that we thoroughly 
understand the systematics in these measurements. This work investigates how the [N ll] fHa ratio based metal- 
licity gradients change with angular resolution, signal-to-noise (S/N), and annular binning parameters. Two ap- 
proaches are used: 1. We downgrade the high angular resolution integral-field data of a gravitationally lensed 
galaxy and re-derive the metallicity gradients at different angular resolution; 2. We simulate high-redshift inte- 
gral field spectroscopy (IPS) observations under different angular resolution and S/N conditions using a local 
galaxy with a known gradient. We find that the measured metallicity gradient changes systematically with 
angular resolution and annular binning. Seeing-limited observations produce significantly flatter gradients than 
higher angular resolution observations. There is a critical angular resolution limit beyond which the measured 
metallicity gradient is substantially different to the intrinsic gradient. This critical angular resolution depends 
on the intrinsic gradient of the galaxy and is < 0.02" for our simulated galaxy. We show that seeing-limited 
high-redshift metallicity gradients are likely to be strongly affected by resolution-driven gradient flattening. 
Annular binning with a small number of annuli produces a more flattened gradient than the intrinsic gradient 
due to weak line smearing. For 3-annuli bins, a minimum S/N of ^ 5 on the [N ll] line is required for the 
faintest annulus to constrain the gradients with meaningful errors. 



Subject headings: galaxies: abundances 
lensing: strong 



galaxies: evolution — galaxies: high-redshift — gravitational 



1. INTRODUCTION 

The chemical abundance distribution within galaxies at lo- 
cal and high redshift offers unique insights into the forma- 
tion and evolution of galaxies. Indeed, the first simple forma- 
tion models of our Milky Way were built upon observational 
knowledge of the Galactic chemical distribution and stellar 
dynamics (Eggen et al. 1962; Searle & Zinn 1978). 

For local galaxies, the existence of radial metal- 
licity gradients has been established in most spi- 
ral galaxies using a bundances in Hii r egions 
jPaael & Edmunds' '198U lEdmunds & Page| 119841; lEvans 
1986; Garnett & Shields 1987; Vilchez et allll988t IShields 
1990; Vila-Costas & Edmundsj |1992[ iZaritskv et all Il994t 
Dutil & Rov 1999). The radial metallicity gradients are 
negative, with a higher metallicity towards the galactic 
center On the other hand, mergers and barred galaxies have 
shallower or flattened gradients due to interaction-induced 
gas flows or bar induced gas inflows (iKewley et al.l 120101: 
iRupke et al. 2010a b). 

In the classic galactic chemical evolution models, the radial 
metallicity gradients formed during inside-out galaxy mass 
assembly. In this scenario, ga s infall timescales i n crease 
with galactocentric distance ("Prantzos & Boissiei^ '20001 
ChiapDini et al. 2001; Molla & Di'az 2005; Magriniet^ 
2007t iFu et al.l l2009). However, even among inside-out for- 



mation models, no consensus has been reached on the direc- 
tion and magnitude of the cosmic time evolution of the metal- 
licity gradient. 
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The time evolution of the metallicity gradient in disk galax- 
ies provides tight constraints on the historical events that mark 
the galactic disc structure evolution. The cosmic metallicity 
gradient evolution is directly linked with the size-growth, ex- 
ternal gas accretion, galactic-scale outflows, and any internal 
transportation of the sta r-forming gas in galaxies. Ear ly an- 
alytical models such as lEdmunds & Greenhowl (119951) have 
provided perceptive predictions on how the gas flows can 
change the origin and evolution of the metallicity gradi- 
ent. Current ACDM cosmological hydrodynamic simulations 
now incorporate the prescriptions for galactic chemical evo- 
lution models and initial predictions for the metallicity gra- 
dient evolution with r edshift are beginning to emerge (e.g., 
iPilkington et alj 120121) . To guide and compare with these 
models and simulations, it is of crucial importance to establish 
an observational baseline for the time evolution of metallicity 
gradients in galaxies. 

Observational constraints on the time evolution of the 
metallicity gradient have been scarce and are tradition- 
ally based on indirect methods using different age-tracers 
such as Cepheid s, planetary nebulae, and open clusters in 
the Milky W ay ( And rievskv etaTI l200l iMaciel etaUlIOOl 
iMac iel & Costa 2010). Direct observations of metalUcity gra- 
dients at high redshift have not been possible until the recent 
employment of high sensitivity near infrared (NIR) integral 
field spectrographs (IFS) on large telescopes. 

As a result, it is only in the past three years that the first ra- 
dial metallicity gradie nt measurem ents have been made for 
high-redshift galaxies dJones et al.li 2010; Cr esci et al.l 1 20101 ; 



Yuan et al] 



2011 



lOuevrel et al.112012; Swinbank et al.l 
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Jones et al]|2012l) . These integral field unit (IFU) studies can 
be divided into three categories according to their angular 
resolution: (1) gravitationally l ensed galaxies with adaptive 
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120111: iJones et all 120121): (2) non-lensed g alaxies with AO 
aided observations (lSwinbanket"ani20T2l): (31 non-lensed 
galaxies with seeing -limited observations TCresci et al.ll2010l: 
lOueyrel et al.ll2012li i. It is interesting to note that the steepest 
gradients discovered so far are all observed in the highest res- 
olution category (1), whereas much shallower gradients are 
reported in categories (2) and (3). 

Compared to local observations, the most obvious restric- 
tions at high redshift are angular resolution and signal-to- 
noise (S/N). Because of these impediments, simplifications 
in data analysis such as annular binning are used to gain S/N 
in spectra. It is unknown whether/how the observational lim- 
itations and data analysis techniques for high redshift stud- 
ies cause systematic uncertainties in the measured metallicity 
gradients. 

With the rapid progress in metallicity gradient observations 
at high redshift, it is imperative that we understand the sys- 
tematics in our measurements before embarking on surveys 
of large samples. 

This paper is the first of a series devoted to understanding 
the systematics of metallicity gradient measurements. Specif- 
ically, this paper investigates how the metallicity gradient is 
affected by the angular resolution, S/N, and annular binning 
that are characteristic for high-z studies. 

We use two approaches: 1. We downgrade the highest an- 
gular resolution and S/N IFU data from our published AO- 
aided lensed galaxy to lower resolutions and then recalculate 
the metallicity gradients at different resolutions; 2. We simu- 
late the IFU data at high-z using a local starburst galaxy with 
a known gradient. We then compare the metallicity gradients 
measured from the simulated data at different angular reso- 
lutions and S/N ratios. We find that the metallicity gradient 
changes systematically with angular resolution, S/N and the 
annular binning parameters. 

This paper is organized as follows. In Section|2] we list the 
assumptions used. In Section [3] and Section HI we describe 
our two approaches and results. We compare our results with 
literature data in Section |5] We discuss our explanations for 
the systematic uncertainties in Section|6] We present our con- 
clusions and future directions in Section |7] Throughout this 
paper we use a standard ACDM cosmology with Ho= 70 km 
Mpc"\ fiAf=0.30, and 17a=0.70. 

2. OVERALL ASSUMPTIONS 

Many factors can affect metallicity gradient measurements. 
In order to analyze the role of angular resolution, S/N and 
annular binning parameter, it is necessary to decouple these 
effects from other factors that may affect the metallicity gra- 
dient. We hold all other parameters that may affect the gra- 
dient constant. To focus on the systematics of the angular 
resolution, S/N and annular binning on the metallicity gradi- 
ent measurement, we find it necessary to impose the following 
assumptions and simplifications: 

• We use the [Nnl/Ha ratio as calibrated in 
iPettini & Pagell (120041) (the PP04N2 method here- 
after) to calculate metallicity. Because of the relatively 
easy access, [N ll] /Ha ratios will continue to be the 
most practical metallicit y diagnostics for hi g h- z st udies 
in the near future. In iKewley & EUisoiil (l2008h . we 
showed that the absolute abundance scale differs for 
different calibrations. However, relative metallicities, 
such as metallicity gradients are robust to within 0.03 
dex on average. To avoid problems with the absolute 



abundance calibration scale, we only apply the single 
PP04N2 diagnostic. We have verified that the results 
for our local galaxy are unchanged when independent 
diagnostics su ch as the [Kobulnicky & Kewley ( 2004) 
and lMcGaughi d 1 99 Ih caUbrations are applied. 

• We assume that the metallicity gradient can be fitted by 
one linear function (c.f., Bresolin et al. 2009b). We ex- 
press gradients in dex per kpc. We do not consider how 
the intrinsic sizes of galaxie s may affect the m etallicity 
gradient measurement (e.g.. lJones et al.ll2012l) . 

• We assume that the metallicity distributions in 
disk galaxies have cylindrical symmetry, i.e., az- 
imuthal abundance variati ons jire not si gnificant (e.g., 
iKennicutt & Gametll996HBresoiin et all2009al) . 

• We assume the galactic center is well-defined and is 
aligned with the peak of Ha emission. Locating the 
galactic center is straightforward for local galaxies, but 
may not be unambiguous for high-z galaxies because of 
their irregular morphology (e.g., Abraham etal. 19961 
iGiavalisco et al. 1996). For similar reasons, we do not 
correct for the inclination angle of the galaxy as it is not 
well constrained for high-z galaxies. Our aim is to sim- 
ulate high-redshift data and to analyze that data under 
the same assumptions and limitations that affect high 
redshift galaxy observations. 

• We only consider metallicity gradients for isolated 
disks. 

• The inhomogeneity of metal distribution has been 
found in both the radial and vertical directions of 
the galactic discs, as well as in early-type galax- 
ies. Although abundance gradients in the vertical 
direction and in early-type galaxies are equally im- 
portant ( Franx & Illingworth 1990: Henry & Worthey] 
ll999i:IMarsakov & Borkova 2006), we confine our dis- 
cussion to the radial gas-phase metallicity in disk galax- 
ies. 

• The uncertainties stemming from l ensing models and 
the effects of shocks and AGN (e.g.,'Wright et al."2010t 
Rich et al. 201 1; Yuan et al. 2012; Westmoquette et ay 
201 2D are beyond the scope of the paper and will be 
investigated in a future paper in this series. 

Note that each of the above items may play an important 
role in metallicity gradient studies. We hold them constant 
in order to filter out the systematics of the angular resolution, 
S/N and annular binning on the metallicity gradient measure- 
ment. By comparing metallicity gradients obtained for the 
same data within the same galaxy using the same calibration, 
we avoid the effects of the issues highlighted above. 

3. METHODOLOGY A: DOWNGRADE HIGH-REDSHIFT HIGH 
ANGULAR RESOLUTION DATA 

The highest intrinsic angular resolution achieved (^ 0.02") 
in metallicity gradient measurements at high-z is in cate- 
gory (1): gravitationally lensed galaxies with adaptive optic s 
(AO) (Jone s etaTI 12010: Yuan etal. 2011; Jo nes et al.ll20T2l) . 
Among these studies, the grand-design spiral Spll49 at z = 
1.49 is most spectacular and suitable for detailed metallic- 
ity analysis because of its clear face-on morphology and for- 
tuitously uniform lensing magnification (Yuan et al.. .201 1,) . 
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Note that all high-z observations of galaxies suffer from the 
loss of low surface brightness pixels. The magnification of 
gravitationally lensed galaxies helps to alleviate the problem 
by bringing more spatial elements with low surface brightness 
within the detection limit compared to non-lensed cases. We 
therefore use Spl 149 as our high redshift testing galaxy. 

3.1. Integral field spectroscopy of Spl 149 
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Fig. 1. — Metallicity ([Nii]/Hq! ratio) vs. radius for Spl 149 downgraded 
to different angular resolution FWHM = 0.1, 0.15, 0.3, 0.75, and 1.2 arcsec 
(in different colors). The solid lines are the best-fit of the metallicity slope, 
the i-cr errors of the slope from linear regression are shown as shaded regions 
around the sohd lines. The imbedded panel on the lower left shows the orig- 
inal OSIRIS Ha map with the three annuli marked. The original data has an 
angular resolution of 0. 1 " (FWHM), corresponding to the pixel scale of the 
data. Note that because of lensing magnification, the corresponding intrinsic 
angular resolution in this case is 0.02". 



Spll49 was observed using the Laser Guide Star Adap- 
tive Optics aided integral field spectrograph (OSIRIS; 
iLarkinet al.ll200 6l) on KECK II in 2011. We achieved an 
image-plane angular resolution of O'.'l, corresponding to a 
source-plane resolution of ^ 0.02" (170 pc). The pixel scale 
we use for Spl 149 is O'.'l, with a field of view (FOV) of 4.8" 
by 6.4" and a spectral resolution of R ^3400. 

Spl 149 is a rare lensing case as it is stretched almost 
equally by ^ 5 times on each side of the two-dimensional 
image. As a result, the source-plane morphology looks al- 
most identical to the image-plane morphology. The error on 
the lensing magnification of Spl 149 is < 10% jSmith et alj 
12009; Yuan etal. 2011.) 

In [Yuan et aLj ( 1201 II) . we derive the metallicity gradi- 
ent for Spl 149 by integrating the spectrum in three an- 
nuli corresponding to a physical length of rl=0.72±0.1 kpc, 
r2=2.34±0.2 kpc, r3=4.5±0.4 kpc. The choice of the annuli 
is such that the integrated spectra could reach S/N > 5 for the 
weak [N ll] lines. The [N ll] fine is robustly detected at S/N > 
5 for the inner two annuli and is a 3a detection for the outer 
annulus which we give an upper limit. In the next subsection, 
we use the same steps to derive metallicity gradients on the 
resolution-downgraded data of Spl 149. 

3.2. Spl 149 Downgraded to Different Angular Resolutions 

For each wavelength slice, we convolve the spatial pix- 
els with a Gaussian kernel of a range of full width half 
maximum (FWHM) values (0.2"-0.8"). We then re-extract 



the spectra using the same three annuli as in lYuan et aTl 
(2011) . We use th e same line-fitting procedures as used 
in .Yuan et al.l (1201 Ih : Gaussian profiles were fitted simul- 
taneously to the three emission lines: [N ll] A6548, 6583 
and Ha. The line profile fitting was conducted using a 
minimization procedure which takes into account the greater 
noise level close to atmospheric OH emissions. The cen- 
ti-oid and velocity width of [N ll] AA6548, 6583 fines were 
constrained by the velocity width of Ha A6563, and the ra- 
tio of [N II] A6548 and [N ll] A6583 is constrained to the the- 
oretical value (Osterbrock 1989). We fit the [Nil] /Ha ra- 
tios in the three annuli and obtain the metallicities using the 
PP04N2 method (Pettini & Pagel2004) . In Figure [T] we show 
the metallicity gradient and 1-a error of the fit to metallicity 
([N ll] /Ha ratio) vs. radius for the downgraded IFU data in a 
few cases of angular resolution FWHM. 

We see from Figure [T| that the measured metallicity gradi- 
ent flattens with poorer angular resolution FWHM. With an 
angular resolution FWHM above ^0.75", inverted (positive) 
gradients begin to appear within the errors of the linear fitting. 

To determine whether the trend continues at even lower an- 
gular resolution, we downgrade the data using a range of an- 
gular resolution FWHM between the original 0. 1 " resolution 
and a fiducially large 2.0" resolution. We find that relation 
between the metallicity gradients and the angular resolution 
FWHM follows an interesting curve (Figure|2l. 




FWHM (orcsec) 



Fig. 2. — The measured metallicity gradient and the l-cr error as a function 
of the FWHM angular resolution based on the data of Spl 149 at z = 1.49 
(Yuan et al. 201 1). The smallest angular resolution con'esponds to the actual 
observed resolution of 0. 1 " or 0.02" intrinsic resolution after correcting for 
lensing magnification. 



Figure |2] suggests that there is a critical angular resolu- 
tion, below which the measured gradient is significantly more 
steepend. Above this critical resolution, e.g., under seeing- 
limited conditions, the gradient approaches a very shallow or 
nearly flat slope. 

The behavior of the gradient vs. angular resolution is wor- 
risome. Ideally, the measured gradient should not be a func- 
tion of angular resolution within the observational errors . If a 
higher angular resolution more closely represents the real gra- 
dient, then what is the angular resolution required to recover 
the intrinsic gradient ? We address this question in Section |4l 

4. METHODOLOGY B: SIMULATE HIGH-REDSHIFT IFU 
OBSERVATION USING LOCAL DATA 
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Fig. 3 . — The continuum and wavelength-collapsed Ho image derived from 
the WiFeS data of Spl7222 (first row), and from the simulated Spl7222 
at 2 = 1.49 in different angular resolution FWHM (row 2-3: examples of 
FWHM = 0.1", 0.2", 0.3"). We see that at FWHM > 0.3", the Ha clumps 
on the outer disk of Spl7222 are not resolved anymore. The artifacts (black 
dots) on the simulated data are regions of low S/N (< 5) on the Ha line. 



4.1. Local Spiral Galaxy IRAS Fl 7222-5953 (Spl 7222 ) 

In order to simulate an IFS observation at high-z, we 
choose a local galaxy that has star formation concentrated at 
the galactic nucleus instead of in the spiral arms. We use 
the isolated luminous infrared galaxy IRAS F17222— 5953 
(Spl7222) aX z = 0.021. Spl7222 is relatively face-on, has a 
morphological type of Sbc, and a moderate SFR (Ha) of ^ 11 
Mq yr^^. The optical IFU data of Spl 7222 is adopted from 
the Wide Field Spectrograph (WiFeS) Integral Field Unit 
(IFU) Great Ob servatory All-Sky LIRG Survey (GOALS) 
Sample (WIGS) (iRich et al.ll20"Tl . 

Briefly, the IFU data were taken using WiFeS at the Mount 
Strom lo and Siding Spri ng Observatory (SSO) 2.3 m tele- 
scope dPopita et al.l2010l) . The blue and red spectra of WiFeS 
were observed with a resolution of R3000 and R7000 and 
wavelength coverage of -3500-5800 A and - 5500-7000 °A 
respe ctively. We use the fully reduced and flux calibrated data 
from IRich et al.1 (l2012h . The pixel scale of WiFeS is 0.5", 
with a field of view (FOV) of 25" x 38". The reduced data are 
binned by 2 pixels in the spatial direction, yielding an effec- 
tive spatial resolution of 1 .0" x 1 .0". 

4.2. Simulating IFU Observations of Spl7222 at z ~ 1.49 



To make realistic comparisons with observed IFU data, we 
fix Spl7222 at the redshift of Spl 149 (z = 1.49). We derive 
the size and surface brightness distribution for each redshifted 
wavelength slice using basic equations of angular diameter 
distance and luminosity distance, and based on the conserva- 
tion of total luminosity. Due to pure cosmology effects, the 
angular size and total flux of Spl7222 are reduced by — 0.05 
and — 2.9 x 10^^ respectively at z — 1.49. 

Because the angular size becomes smaller at higher red- 
shift, an angular scale of 0.02" is required to fully recover 
the spatial samplings of the local IFU data. We thus set the 
highest FWHM resolution of our simulation as 0.02". 

The original Poisson-noise dominated optical spectra are 
shifted into the NIR at z = 1.49, and real NIR data are sky- 
background dominated. We compose a noise spectrum at the 
spectral resolution of WiFes by interpolating the sky -residuals 
of our OSIRIS datacube. The noise spectrum represents the 
sky OH emission residuals in a real observation in the NIR 
at z = 1.49. The noise spectrum is scaled and added to the 
redshifted spectrum of each spaxel according to the required 
input S/N. 

We define the input S/N of our simulation on the spaxel 
where Ha emission peaks. For example, an input S/N of 100 
means that the noise spectrum is scaled and added to the target 
spectrum such that the S/N ratio of the Ha line on the bright- 
est pixel is 100. To facilitate comparison with observations, 
we also measure the observed S/N of the [N n] lines in outer 
annulus on the simulated datacube (Section 4.5). 

Note that in this simulation, we do not consider any intrinsic 
evolutionary effect except that we multiply the total flux of the 
original IFU data by a factor of 20. The factor of 20 is adopted 
for two reasons: 1. The SFR of galaxies at z 2 are found 
to be — 20 higher than galaxies at z ~ at a fixed mass o f 
- 10^° Mq (e.g.. Noesk e et^HMfM iZahid et al.ll27)T2h : 
2. Without manually increasing the total flux, Spl7222 would 
have been under the detection limit of any current NIR IFU 
instrument. We find that by increasing the intrinsic flux by 20, 
Spl 7222 would be well detected in NIR IFU spectrographs 
such as OSIRIS on KECK. The factor of — 20 also coincides 
with the lensing flux magnification of Spl 149 in Section[3] 

Finally, following similar steps as in Section[3] we convolve 
the spatial pixels with different angular resolution FWHM and 
generate a set of simulated IFU data at z = 1.49. 

For both the local and simulated data cube, we use the 
same line-fitting procedures as described in lYuan et al.l(l201 1[) 
and in Section 3.2. Figure |3] shows the continuum and Ha 
emission line image derived from the original WiFeS data of 
Spl7222, and from the simulated high-z Spl7222 datacube in 
different angular resolutions. 

4.3. The Effect of Annular Binning on Metallicity Gradient 
Annular binning is a commonly used method to derive 



metallicity gradients fo r low S/N high-z data (Jone s et al 
| 2010- Yuan et al. 20111: lOueyrel et al. 2012: Swinbank et al 
I2OI2I) . In order to quantify the effect of annular binning on 
metallicity gradient measurements, we compare the effect of 
3-annuli binning with full sampling without binning. Since 
only the local datacube has sufficient spatial resolution and 
S/N simultaneously to compare a full and 3-annuli sampling, 
we test these 2 types of samplings on the original datacube of 
Spl 7222 first, without degrading the resolution. 

We show the metallicity gradients derived from a 3-annuli 
binning and a full sampling in radius in the lower right panel 
of Figure |4] The three annuli are defined on the simulated 
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Fig. 4. — Derived [N 11] /Ha -based metallicity gradient for the local data of Spl7222. Left: Ha intensity derived from the WiFeS datacube. We use the same 
three annuli as defined on the simulated high-z data. The average spectra (zoomed in for the [N 11] and Ha vicinity) in the three annuli are shown on the light 
panel (al-a3). The flux is in unit of 10~^^ ergs s~^ cm~^ A^^ . The panel on the bottom right shows the metallicity gradient derived using the three annuli 
(black) in comparison with a full sampling (red) without binning. 
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Fig. 5. — The percentage of the difference between the measured metallic- 
ity gradient in annuli and the true metallicity gradient in full sampling as a 
function of the number of annuius used. The blue circles are the measure- 
ments from the annuli that are unifonnly divided from the galactic center to 
the outer disk. The red circle shows the measurement from the 3 annuli that 
are defined on the high-z data constrained by the S/N of the weak [Nil] lines 
(Section 4.3). 



high-z data and they are chosen such that the [N n] line from 
the outer-most annuius is robustly detected (S/N > 5) on all 
annuli of the hig h-z datacube, consistent with current high-z 
studies jJones et al. 2010; Yuan et al. 2011). 

We see from Figure|4]that the gradient from annular binning 
is shallower than the gradient derived from full samplings. In 
the case of Spl7222, binning in 3 annuli yields a metallic- 
ity gradient of -0.016±0.004 dex kpc^^, whereas the fulling 
sampling yields -0.032±0.003 dexkpc^^. 

To further investigate the relation between the measured 
gradient and the number of annuli used, we calculate the mea- 



sured metallicity on the local IFU data of Sp 17222 using a 
range of annuli. Note that the annuli in Figure |4] are defined 
on the simulated high-z data and are constrained by the S/N 
of the [N II] line from the outer-most annuius. The choice of 
the annuli is therefore defined on the emission line flux dis- 
tribution of the high-z data. A more straight-forward method 
of defining annuli is to divide the radius equally into N uni- 
form bins. These uniform annuli can be applied to the lo- 
cal data where the requirement of S/N> 5 on the [Nil] line 
is not a constraint for most radii. We therefore divide the 
Spl7222 data equally in radius using a number of N,„i„ = 3 
and N„iax = 20 annuli. We show in Figure|5]the percentage of 
the difference between the measured metallicity gradient from 
different numbers of annuli and the true metallicity gradient 
from full sampling as a function of the number of annuli. We 
find that with N=3-5 uniform annuli, the measured gradients 
are close to the true gradient within ~ 30%; with N>6 uni- 
form annuli, the true gradient can be recovered to within 10% . 
We also show that the metallicity gradient from the uniform 
3-annuli is closer to the true gradient (28% c.f. 54%) than the 
S/N constrained 3-annuli employed on the high-z data. 

In summary, we see that binning in 3 (or another small num- 
ber of) annuli introduces non-negligible systematic errors on 
the metallicity gradient measurement. The definition of the 
location of the annuli also plays a non-negligible role in the 
metallicity gradient measurement. The choice of the annuli 
used in high-z studies is weighted by the S/N of the weak 
[N II] lines and can only recover the true gradient to within 
54% in the simulated case of Spl7222. The choice of uni- 
formly distributed annuli recovers the true gradient signifi- 
cantly better In the case of Sp 17222, using N > 6 annuli 
can recover the true gradient to within 10%. 



4.4. The Effect of Angular Resolution on Metallicity 
Gradient Measurements 
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Fig. 6. — Examples of the [N II] and Ha lines exti'acted from the 3 annuli on a simulated datacube (S/N=100, FWHM=0.02"). The panels are organized in the 
same order as in Figure |4] with the same flux and wavelength units. The panel on the bottom right shows the metallicity gradient (blue) based on the [Nil] and 
Ha line ratios. The local 3-annuli (black) and full sampling (red) gradients are also shown for comparison. 
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Fig. 7. — Left: Metallicity vs. radius on the simulated datacube (S/N fixed at 100). Three-annuli bins are used. We show examples of a few angular resolution 
FWHM in different colors. For compaiison, the three annuli and fully sampled metalhcity of the local data are shown in red. Right: Simulated metallicity 
gradient as a function of angular resolution FWHM. The horizontal red lines show the location of the metallicity gradients measured on the local data using three 
annuli and full samplings respectively. We see that the local gradient is not recovered even in the highest resolution of the simulation (FWHM=0.02 "). 



We extract the 3-annuli spectra for our simulated high-z 
datacube with a range of angular resolution FWHM and S/N 
ratios. Figure |6] shows an example of the 3-annuli spectra and 
derived metallicity gradient. The example in Figure |6] has an 
angular resolution of FWHM=0.02" and a limiting S/N of 100 
on the Ha line. A S/N of 100 on the Ha line corresponds to a 
measured S/N of 7.5 on [N n] line of the faintest annulus. 

Figure |7] shows the PP04N2 -based metallicity gradient for 
different angular resolution FWHM. We see the same trend 



in the simulated data as seen in Figure [T] the measured 
metallicity gradient is shallower as the FWHM increases. As 
the angular resolution FWHM approaches the seeing limited 
regime, the gradients are essentially consistent with a flat 
0) value. The curve of the metallicity gradients vs. the angu- 
lar resolution FWHM relation is similar to Figure |2] However, 
the change of the FWHM from steep to shallow gradients is 
slower. 

The deviation of the measured gradient from the true gra- 
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Fig. 8. — The deviation of tlie measured gradient from tlie "true gradient" as 
a function of S/N (angular resolution fixed at 0.02"). The horizontal red line 
shows the true gradient measured on the fully sampled local data of Spl7222. 
The S/N marked on the bottom are defined on the peak spaxel of the Ha 
emission line (Section 4.2). The corresponding S/N of the [N II] line on the 
faintest annulus are marked on the top. 
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Fig. 9. — The fact that measurements for different galaxies in different 
studies follow the same trend as Figure 2 and 7 strongly suggest that the 
angular resolution has played a role in causing the systematic difference in the 
metallicity gradient. The effect of angular resolution needs to be quantified 
before interpreting the metallicity gradient evolution with redshift. 

dient measured from the fully sampled local data of Sp 17222 
can be derived from Figure |7] The closest gradient measure- 
ment we can obtain is 0.095 ± 0.004 dex perkpc (true gradient 
= 0.032±0.004 dex per kpc), i.e., only 29% of the true gra- 
dient is recovered by the highest angular resolution and S/N 
observation simulated. 

4.5. The Effect of S/N on Metallicity Gradient Measurements 

Section 133] shows the high input S/N (100) case. As de- 
scribed in Section|42] the input S/N is defined on the Ha line 
of the central pixel. 

We repeat our analysis in Section |4~4l using a range of i nput 
S/N. We vary the S/N ratios as described in Section 14.21 We 
fix the angular resolution at FWHM=0.02". Figure |8] shows 
the deviation of the measured gradient from the true gradi- 
ent as function of S/N. The errors of the metallicity gradient 
are derived from the errors in the slope of the two-variable 



(metallicity and radius) linear regression. We see that the most 
significant effect of S/N is on the errors of the metallicity gra- 
dient measurement. The errors of the gradients decrease with 
S/N. At Ha S/N (peak) < 50 or [N ll] S/N (3'''^ annulus) <5, 
the error bars of the gradient are so large that the gradients 
can not be meaningfully constrained, i.e., both a positive and 
negative gradient are consistent within the errors. 

5. METALLICITY GRADIENT VS. ANGULAR RESOLUTION FROM 
THE LITERATURE 

As described in Section 1, current IFS studies on metallic- 
ity gradients can be divided into three categories according 
to their angular resolution: (1) gravitationally lensed galax- 
ies with AO; (2) non-lensed galaxies with AO; (3) non-lensed 
galaxies under seeing-limited conditions. We gather the lit- 
erature data in these 3 categories and examine the measured 
metallicity gradients with respect to the angular resolution 
FWHM used in the observations. 

For category (1 ) we use t h e 5 l e nsed galaxi e s from 
iJones et al.1 (12010) ; lYuan et al.1 (l20Tll) ; iJones et alJ JMl) . 
We exclude 1 merging system from the lensed galaxies be- 
cause galaxy mergers are known to flatten metallicity gra- 
dients dKewlev et alj [loTol; iRupke et alj 1201 Obi; iRich et alj 
[Ml. For (2} . we u se the 7 AO-correcte d galaxies from 
'Swinb anket alJ (120121) . The galaxy sample in lSwinbank et alJ 
(2012) is chosen from the Hig h-Z Emission Line Survey 
(HiZELS) survey (ISobralll2012b which target s Ha emitters 
close to bright stars. For (3), we use the data of | Quevrel et"ar] 
(120121) . as 17 out of the 18 galaxies in lOuevrel et al.l (l2012h 
are seeing-limited measurements, with 1 AO corrected mea- 
surement. The galaxy sample in Ouevrel et al. (2012) is 
part of the MASSIV project ( Contini et al.ll2012l) and is cho- 
sen based on the visibility of the Ha line in J or H band. 
These 1 8 galaxies do not inc lude interacting galaxies iden- 
tified in 'Ouevr el etall ( 120121) . Since [N ll] and Ha Unes 
are available for all the samples, we calculate the metallici- 
ties using the PP04N2 metallicity diagnostic. Note that the 
metallicity gradients in jQ uevrel et al. (2012) are based on the 
iPerez-Montero & Continii (2009) [Nil] /Ha metallicity cali- 
bration, which we have converted to the PP04N2 metallicity 
calibration. The metallicity gradients are shown as a function 
of the angular resolution FWHM in Figure |9] 

Note that there is a redshift difference among the three sam- 
ples. The median redshift for sample (1) is z ^ 2, the median 
redshifts for samples (2) and (3) are z ^ 1.5 and z ~ 1.3 
respectively. If the angular resolution effect does not play a 
role in these observations, one would be tempted conclude 
that the metallicity gradient is evolving from the steeper gra- 
dients at z ^ 2 to much shallower gradients at z ^ 1. How- 
ever, relating Figure |9]with Fi gurel2]and Figure |7] we see a 
striking similarity in the behavior of the measured metallicity 
gradient vs. the angular resolution FWHM relation. The fact 
that the steepest gradients are only seen in the lensing-i-AO 
samples (including one galaxy at 1.5) is consistent with 
our findings in the previous sections that there exists a criti- 
cal angular resolution limit beyond which the absolute value 
of the gradients are significantly under-estimated. It is possi- 
ble that angular resolution causes the trend in Figure |9] Until 
the effect of angular resolution on metallicity gradients is well 
understood, extreme caution must be taken when interpreting 
apparent evolution in metallicity gradients with redshift. 

6. DISCUSSION 
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6.1. Explanation of the angular resolution effect: "low S/N 
line smearing" effect 

The PP04N2 based metallicity depends on the ratio of the 
[N II] and Ha lines. The [N ll] lines are usually more than 3 
times weaker than Ha lines for normal star-forming regions. 
Given a negative radial metallicity gradient, much weaker 
[N II] lines are expected in the outer regions of a galaxy. Any 
spatial averaging/smoothing process selectively smears the 
low S/N [N II] line regions. As a result, the spatially averaged 
spectra are weighted more towards the regions of stronger 
[N II] lines, leading to over-estimation of the metallicity in 
the outer-disk. Accordingly, binning using small numbers of 
annuli has the same effect as lower resolution (i.e. larger an- 
gular resolution FWHM). 

Using this argument, we would expect that the steeper the 
intrinsic gradient is, the larger the observed gradient is going 
to deviate from the intrinsic value. The curve of the gradient 
vs. FWHM relation (e.g.. Figure |2] Figure |7] and Figure |9) 
would depend on the intrinsic shape of the metallicity gradi- 
ent. 

If low S/N line smearing is the cause of the angular resolu- 
tion effect observed in Figure |2] Figure |7] and Figure |9] one 
would think that it is possible to cancel out the smearing ef- 
fect by using a reversely weighted function. However, it is 
not straightforward to establish this anti-weighting function, 
as it relies on knowing the intrinsic gradient a priori. More- 
over, if the "smearing" process is done by the atmosphere, 
it is impossible to de-convolve the observed FWHM back to 
the pre-smeared values, as it would be equivalent to finding 
an algorithm to "recover" the data with any wished angular 
resolution. 

7. CONCLUSION 

In this paper, we have demonstrated that the angular resolu- 
tion, S/N and annular binning introduce significant systematic 
errors in metallicity gradient measurements at high-z. 



We find that; 

• The measured metallicity gradient changes systemati- 
cally with angular resolution FWHM. Seeing-limited 
observations are likely to produce more flattened gra- 
dients than AO-aided high-resolution studies. 

• There is a critical angular resolution FWHM range be- 
yond which the measured metallicity is significantly 
more flattened than the intrinsic metallicity. This crit- 
ical FWHM depends on the intrinsic gradient of the 
galaxy. For the two cases used in this work, the crit- 
ical FWHM is < 0.02", only currently reachable with 
AO + gravitational lensing. 

• For a fixed angular resolution, the errors of the metal- 
licity gradient increase as S/N decrease. At low S/N 
(< 5 for [N II] line in the simulated case), the errors 
are so large that the gradient can not be meaningfully 
constrained. 

• Three-annuli binning or any limited number of annular 
binning yields a more flattened gradient than the intrin- 
sic gradient. 

Until these effects are thoroughly understood, we urge cau- 
tion in interpreting metallicity gradient evolution with red- 
shift. Our next work is to build an empirical angular reso- 
lution library by simulating a large sample of high-z galaxies 
using local galaxies with known gradients. 
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